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Q . We study the progenitor dependence of the black hole formation and its as- 

sociated neutrino signals from the gravitational collapse of non-rotating mas- 
^ ! sive stars, following t he preceding study on the single progenitor model in 



Sumiyoshi et al. ( 200? !). We aim to clarify whether the dynamical evolution 



toward the black hole formation occurs in the same manner for different pro- 
genitors and to examine whether the characteristic of neutrino bursts is general 
having the short duration and the rapidly increasing average energies. We per- 
form the numerical simulations by general relativistic z/-radiation hydrodynamics 
to follow the dynamical evolution from the collapse of pre-supernova models of 
4OM0 and 5OM0 toward the black hole formation via contracting proto-neutron 
stars. For the three progenitor models studied in this paper, we found that the 
black hole formation occurs in ^0.4-1.5 s after core bounce through the increase 
of proto-neutron star mass together with the short and energetic neutrino burst. 
We found that density profile of progenitor is important to determine the ac- 
cretion rate onto the proto-neutron star and, therefore, the duration of neutrino 
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burst. We compare the neutrino bursts of black hole forming events from differ- 
ent progenitors and discuss whether we can probe clearly the progenitor and/or 
the dense matter. 

Subject headings: supernovae: general — stars: neutron — black hole physics — 
neutrinos — hydrodynamics — equation of state 



1. Introduction 



The final fate of massive stars is one of key issues in stellar physics as well as astrophysics. 
The massive stars beyond ^SMq are the origin of compact objects, i.e. neutron stars and 
black holes, as the outcome of gravitational collapse at the end of their lives. They contribute 
also to the production of heavy elements through the explosive nucleosynthesis and to the 
high energy phenomena with the bursts of electro-magnetic radiations and neutrinos, which 
in turn play important roles in the evolution of galaxies. The massive stars in the wide 
range between '^SM© and ^IOOMq lead to various explosive phenomena s uch as core-collaps e 
supernovae, hypernovae (or gamma ray bursts) and failed supernovae (IHeger et al.ll2003l ). 
Such energetic phenomena are dependent on the properties of progenitors such as masses, 
rotations and profiles, for example, as in the mass limit of ordi nary supernovae and th e 
branches of hypernovae and failed supernovae as a function of mass ( IMaeda fc Nomotdl2003l ) . 
We focus here on the fate of non-rotating massive stars, as an extreme of failed supernovae, 
in the mass range of ^4O^5OM0. 

Such a death of massive stars attracts attention recently since it is one of possible 
channels of the black ho l e formation and is characte rized by the unique neutrino signals 
(ILiebendorfer et al.l l2004l : ISumiyoshi et al.l l2006l . 120071 ) : the short duration and increase of 
average energies and luminosities in time. This is in contrast to the ordinary supernova 
neutrinos, whose emissions from a nascent pro to-neutron sta r last for ^20 s with a gradual 
decrease of average energies and luminosities (j Burrowsl 1 19881 : ISuzukil 1 19941 ) . The black hole 
formation will be inevitable for massive stars beyond ^4OM0, which was estimated with 



large uncertainties to be the mass threshold for black hole formation by iFryerl (119991 ). The 
dynamical collapse to a black hole is triggered by the intense accretion of matter onto the 
proto-neutron star. Note that this is different from the so-called delay e d black hole for- 
mation of the proto-neu tron star with a fixed mass (IKeil fc Jankal 1 19951 : IPons et al.l Il999l : 
Baumgarte et al.l Il996al ) owing to the termination of accretion after a successful launch 
of shock wave. It should be also emphasized th at the neutrino signal may be used as a 
probe into the equatio n of state of nuclear matter (ISumiyoshi et al.ll2006l ) and exotic matter 
(INakazato et al.l 12008 ) in a different way from the delayed collapse of proto-neutron stars 
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2001 
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This article is the second paper, foUowing the first paper (ISumiyoshi et al.ll2007l ). of serial 
studies on the black hole formation associated with the short burst of energetic neutrinos. In 
the first paper, the evolution toward the black hole formation from the gravitational collapse 
of the 4OM0 star have been examined by taking account of the dependence of equation 
of state. By solving general relativistic ^/-radiation hydrodynamics, they have followed the 
hydrodynamics up to the moment of black hole formation and the neutrino distributions 
during the whole evolution to obtain the detailed features of neutrino emission. They have 
shown that the black hole formation occurs in ^1 s after the core bounce via proto-neutron 
star evolution, but it depends crucially on the equation of state of dense matter. The 
equation of state (EOS) controls the timing of the re-collapse from the proto-neutron star 
to the black hole since the EOS determines the maximum mass. The mass of proto-neutron 
star increases due to the accretion of material in the situation of failed shock wave and the 
dynamical collapse is triggere d when the mass rea c hes th e critical mass. We remark that the 
progenitor model of 4OM0 by Woosley fc Weaver dlQQsh has been unique ly adopted in their 
studies and other studies on such massive stars ( ILiebendorfer et al.ll2004l ). 



In the first paper (ISumiyoshi et al.l 120071 ). associated neutrino burst has been shown to 
be unique to identify the black hole formation and to probe the properties of dense mat- 
ter. The neutrino burst is terminated by the black hole formation and the average energies 
of neutrinos increase toward the end because of high temperature realized in the collaps- 
ing proto-neutron stars. These features are again sensitive to the properties of equation of 
state, which determines the duration of neutrino burst and the properties of compressed 
matter at neutrino emitting region. The neutrino burst toward the black hol e formation 



can, in principle, be detectable at the terrestrial neutrino detector facilities (llkeda et al 



20071) and can be used as a signature of black hole formation and a probe of dense matter 
( ISumiyoshi et al.ll2006l ). Moreover, if a collapse of massive star occurs in a nearby galaxy, 
the e vent may be observed, for example, in the proposed survey of disappearance of massive 
stars (IKochanek et al.ll2008l ) and it will be possible to obtain the information on the progen- 
itor. However, one has to provide the detailed features of neutrino burst by taking care of 
uncertainties of progenitor models in addition to dense matter if one wants to use them as 
templates for the neutrino search. 

In this study, we make a first attempt to see wh ether the features of black hole formation 
reported in the first paper (ISumiyoshi et al.l 120071 ) are taken over to different progenitors. 
We adopt three stellar models as the initial condition of the numerical simulations to make 
the comparison with the evolut i ons b y the numerical modeling we have adopted in the 
previous study (ISumiyoshi et al.l 120071 ). In addition to the stellar model of 4OM0 star by 
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Wooslev fc Weaverl (Il995r >. we adopt the two stellar models bv iHashimotd (119951) and by 
Tominag a et all J2007h . As a case of different mass, we adopt a 5OM0 star bv lTominaga et al.l 



(j2007 ). who studied the evolution of metal-poor massive stars. We have chosen the model 
with zero metallicity among their series of different metallicities. This model is the limiting 
case of metal-poor massive stars without mass loss. As a model for the same mass (4OM0) 
but by a d i fferen t method of stellar evolutionary calculation, we adopt a 4OM0 star by 
Hashimoto (j 199 51). The comparison between the different models enables us to explore the 



dependence of numerical results on the uncertainty in the stellar evolution. 

The investigation on the progenitor dependence is important for the utilization of neu- 
trino bursts as a probe of equation of state. If the neutrino burst is quite different depending 
on the progenitors, it becomes difficult to distinguish the differences found by the equation 
of state. Using the two sets of EOS (LS-EOS and Shen-EOS), it has been shown that the 
time duration of neutrino bursts differ by a factor of ^2 depending on the softness of EOS. If 
the progenitor models give large differences, they may smear out the EOS difference, which 
we want to use as a probe. 

Our aim of the paper is, therefore, to make a first exploration to the dependence on the 
progenitors by adopting a small number of stellar models and to show that the dynamics 
toward the black hole formation is common in the massive stars having a large central 
iron core. We aim also to demonstrate that the neutrino signal has similar features once 
we set the equation of state. We analyze how the evolution is related with the density 
profile of progenitor, which determines the accretion rate. We discuss the characteristics 
of neutrino bursts obtained by different progenitors with the sets of equation of state for 
possible detection in future. 

We plan this paper as follows. We briefly explain the method of numerical simulations 
in ^ We compare the progenitors adopted for initial models and discuss the differences 
in §31 We report the numerical results on the dynamics in § §4.11 and 14.21 for two models 
through the comparison with the previous result. We describe the accretion rates and the 
profiles during the evolution in § §4.31 and 14.41 to explain the differences among the models. 
We discuss the characteristics of neutrino signals in the models to distinguish the progenitor 
and the equation of state in §4.51 Implications of this work and summary will be made in 
§g5]and[6l 
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Numerical Simulations 



Numerical simulations are performed by the nume rical code of general relativistic v- 



radiation hydrodynam ics under the spherical symmetry (jYamadalll997l : lYamada et al 



1999: 



Sumiyoshi et al.l 120051 ). This numerical code has been applied to the gr avitational collapse 
of massive stars to study the supernova explosion (ISumivoshi et al.ll2005l ) and the black hole 
formation dSumivoshi et alJl20od . [2007l : l^ et al.lbooj ^ 120071 ) . The detailed description 

on the numeric al simulations of the gr avitational collapse toward the black hole formation 
can be found i nlSumivoshi et al.l (120071). which reported also the results for the progenitor of 
a 4OM0 star by Woosley fc Weaver! ( 1995 ). We adopt the same number of grids for lagrangian 
mass coordinate, variables in neutrino distributions for z/g, z/g, z/^/^ and P^/^. The mesh for 
mass coordinate is optimized for each progenitor model to follow the accretion phase with 
an enough resolution using the rezoning method. 

The fully implicit method of z/-radiation hydrodynamics enables us to follow the long 
evolution over '^l s toward the black hole formation. The exact treatment of z/-transfer to- 
gether with hydrodynamics in general relativity permits the detailed evaluations of the neu- 
trino fluxes from the contracting proto-neutron stars until the black hole for mation. We de- 



fine th e bla ck hole formation by fin ding the apparent horizon as explained in iNakazato et al. 
(I2OO6I ) and ISumivoshi et al.l (l2007h . 



We adopt the microphysics in the same manner as in lSumiyoshi et al.l (120071 ) to facilitate 
the comparison of new results with the previous result for 4OM0. We use the two sets of 
equation of state (EOS) fo r supernova simula tions to assess the infiuence of EOS on different 
progenitors. The EOS by Shen et al. ( 1998al lbl) (Shen-EOS), which is obtained by the rela- 
tivistic mean field framework with the data of unstable nuclei, provide s a rather stiff EOS 
having the maximum mass of neutron stars to be 2.2M0. The EOS by iLattimer fc Swesty 
(jl99ll ) (LS-EOS) is obtained by the extension of the compressible liquid model, taking a 
density dependence of energy based on the Skyrme-interaction. We have chosen the case of 
incompressibility of 180 MeV among the three choices in LS-EOS as a representative of soft 
EOS having th e maximum mass of neu tron stars to be l.SM©. The same set of weak interac- 



tion rates as in 



Sumiyoshi et al.l (l2007l) i s adopted for the current simulation s by implement- 



ing t he standard formulation of ISruennl (119851 ) and the several extensions (ISumiyoshi et al 
2o3). 
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3. Initial Models 

We adopt three different profiles from the pre-supernova models of massive stars as initial 
models. We have picked massive stars around 4OM0, which are massive enough to have the 
black hole formation without the successful explosion. The choice of the progenitor models 
in the current study is meant for a first and small-scale trial exploration and does not refiect 
a possible great variety of massive stars. We use the central part of pre-supernova profiles 
up to 3M0 in baryon mass coordinate and set density, electron fraction and temperature 
distributions as initial configurations for numerical simulations. 



As a base line model, we adopt a 4OM0 star from the set of massive stars by lWoosley fc Weaver 



( I1995I ). This model has been used for the previous calculation s of our studies o f black hole 
formation and the details of numerical results can be found in lSumiyoshi et al.l (12003). We 
name this model W40 and will use as a reference to examine the new results. The size 
of iron core in W40 model is I.OSMq. As an example of other types of massive stars, we 
adopt a case of 50Mr7^ with zero metallic i ty (Z=0) from the se ries of massive stars having 
low metallicities ( lUmeda fc NomotQll2005l : iTominaga et al.ll2007l ). These zero-metal stars do 
not ex perience the mass l oss during the stellar evolution and remain massive as in the initial 
stage (jUmeda et al.ll2000 l) . This choice is meant to explore the dependence of the phenomena 
on the progenitor mass and the stellar evolution. The size of iron core is I.SSMqUI We name 
this model T50. In order to examine the difference due to the modeli ng of stellar evolutions 
for the same stellar mass, we adopt a 4OM0 star by iHashimotd (119951 ). The size of iron core 
in this model is I.88M0, which is largest among the series of his models up to 7OM0. We 
name this model H40. 

We compare the profiles of the three adopted models in Fig. [H We can see quantitative 
differences among the three models though the profiles look similar at first glance. The 
density of T50 and H40 are higher than those of W40 in the central region whereas the 
order is reversed at the outer region. This difference in gradients of density profile from the 
center to the outer part is related to the difference in accretion rates onto proto-neutron 
star. Actually, the densities at M5 = 2-3M0, which cover the range for the maximum mass 
of proto-neutron stars, are crucial to determine the final stage of accreting proto-neutron 
star through the free-fall time scale \j\fGp). 

The profile of electron fraction in T50 is slightly different from those in W40 and H40. 
There are several kinks at different positions due to the shell structure. The temperature of 



^This value is determined by the criteria of Fg < 0.50, which corresponds roughly to the ^^Ni abundant 
region. A smaller value, 1.46M0, is obtained by the criteria of < 0.49, which is the convention in their 
study (N. Tominaga, private communication, 2006). 
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T50 and H40 are higher than that of W40 in the central region. This difference, however, 
smears out during the coUapse when we compare the profiles at the timing with the same 
central density. We will discuss the profiles later in §4.41 to explain the behavior of accretion 
rates. 



4. Numerical Results 

We present the numerical results for the cases of T50 ( §4.11) and H40 ( §4.21) through the 
comparison with the base line model of W40. We add a character, S or L for Shen-EOS or 
LS-EOS, respectively, in the model names to denote the choice of the EOS. We summarize 
the calculated models in Table [H We discuss the differences among the three models and 
their implications in § §4.3I - I¥y5l 

4.1. SOMq star by Tominaga et al. 

We show the radial trajectories of mass elements as a function of time after bounce 
for models with T50S and T50L in Figs. [2]and[3l respectively. The trajectories are plotted 
for each O.O2M0 in mass coordinate. Thick lines denote the trajectories for O.5M0, I.OMq, 
I.5M0, 2.OM0 and 2.5M0. Thick dashed line represents the position of shock wave. 

The general feature of evolution is found to similar to the case of corresponding W40 
models. After the core bounce, the shock wave is launched up to ^150 km in 0.1 s and 
then recedes gradually toward the surface of central object. The nascent proto-neutron star 
shrinks due to the mass increase slowly in T50S (quickly in T50L) by the intense accretion. 
The dynamical collapse occurs when the proto-neutron star mass exceeds a critical mass, 
which depends on the stiffness of EOS, and the black hole is formed. 

The duration from the core bounce to the black hole formation is 1.51 s and 0.51 s for 
T50S and T50L, respectively, refiecting the difference between SH-EOS and LS-EOS. The 
duration in T50S is slightly longer than that in W40S, while the duration in T50L is slightly 
shorter than that in W40L. This is because the accretion rate is different between T50 and 
W40 in a time-dependent manner as we will see in §4.31 

The duration of neutrino bursts is determined by the time toward the black hole for- 
mation. In Figs. [4]and[5l we display the time evolution of average energies and luminosities 
of neutrinos emitted during the evolution described above. The corresponding results with 
the same EOS but using the progenitor of W40 are shown by thin lines in order to see the 
progenitor dependence. In general, the features of neutrino emission for T50 resemble those 
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for W40. The average energies for all species show the increasing tendency after bounce. 
The persistent increase for v^j^ is remarkable among them. Looking more closely, small 
differences appear at late phase (after ^^0.4 s). The luminosities for T50S become clearly 
lower than those for W40S at the late phase. This is related with a lower accretion rate in 
T50 than that in W40 at ~1 s after bounce (see Fig. [8]). 

We remark that a rapid decrease of neutrino energies by the gravitational redshift that 
is expected at the terminal phase has not been observed yet in the current simulations. One 
has to follow the evolution further (^20 ms) after the for mation of the apparent horizon. 



For that purpose, however, a singularity avoiding scheme (iBaumgarte et al.lll996bl ) should 
be implemented. We note also that some wiggles and/or jumps seen in Figs. H] and [5] are 
numerical artifacts caused by insufficient numerical resolutions at the final phase of numerical 
simulations. 



4.2. 4OM0 star by Hashimoto 

We display the radial trajectories of mass elements as a function of time after bounce for 
model with H40L in Fig. [61 The notation is the same as in Fig. [21 In this case, the evolution 
is much quicker than the cases in W40L and T50L. The time till the black hole formation is 
0.36 s in H40L. Since the accretion rate in H40 is larger than those in W40 and T50, the mass 
of proto-neutron star reaches the critical mass earlier (see §4.31) . The contraction of central 
object is quick accordingly and, therefore, the density and temperature inside become high 
early on. It is recognizable that the position of shock wave behaves in a different manner 
from the other cases. The shock wave stays around 100 km after the initial launch instead 
of the recession. This is mainly because the neutrino luminosity is high due to the high 
accretion rate. 

We show the time evolution of average energies and luminosities of neutrinos in Fig. [71 
The duration of neutrino burst is clearly shorter than that of W40L. The average energy of 
Vfj^/r is larger and increases much faster than that of W40L, refiecting the different temper- 
atures, whereas those of and behave similarly. The neutrino luminosities of all species 
in H40L are higher than those in W40L. The luminosities of and are proportional to 
the accretion luminosity (^ GMM/r) and are high due to a high accretion rate (M) and 
a more massive object (M) at center. The high luminosity of z/^/^ originates from a high 
temperature in the massive object. The fact that the luminosity refiects the behavior of the 
accretion luminosity was also s een for the neutrino em issions at the shock breakout in the 



core-collapse of 11-2OM0 stars (IThompson et al.ll2003l ). The larger accretion rates found in 



the current study of 4O-5OM0 stars lead to more rapid increase of luminosity than for less 
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massive stars. 



4.3. Accretion rate 

As we have discussed already, the accretion rate is a key factor to determine the behavior 
toward the black hole formation. In other words, the increasing rate of proto-neutron star 
mass is crucial to determine the growth time to the critical mass. We show in Fig. [8] the 
baryon mass of central object (corresponding to the proto-neutron star) as a function of 
time after bounce. We define here the central object, which is quasi- hydrostatic, as the one 
inside the position of shock wave. Note that the time derivative of the mass curve (i.e. the 
gradient of the curve) is the accretion rate. After the core bounce, the central object acquires 
^1.5M0 within 50 ms. The proto-neutron star with a typical mass is formed immediately 
after the core bounce. This feature is common for all models of black hole forming collapse 
and is different from the case in or dinary supernovae, where t he proto-neutron star (^l.SM©) 
formation takes place in ^0.3 s (i Burrowsl 1988 : Suzuki 1994 ). The mass increases gradually 



toward the critical mass in each model and the dynamical collapse to black hole occurs at 
the end point. The end point depends on the adopted progenitor and the equation of state. 

Amazingly, the increase of mass in the cases of W40 and T50 are similar each other. 
Although the masses for T50 are slightly larger up to 0.6 s after bounce, their gradients (i.e. 
accretion rate) are the same, judging from the parallel curves between T50S and W40S (T50L 
and W40L). This similarity is the reason why the evolution of T50 and W40 and resulting 
neutrino signals are quite similar. Adopting LS-EOS, the black hole formation occurs around 
the same time (^0.5 s) in the cases of T50L and W40L. In the case of SH-EOS, however, 
there is a difference of about 0.2 s between the cases of T50S and W40S. This is due to a 
smaller accretion rate in T50 at late phase (^0.6 s after bounce) than that in W40 as seen 
in the crossing of curves for T50S and W40S. We note that the small accretion rate at the 
late phase is related with the density profile of outer layer in the adopted progenitor. 

It is remarkable that the behavior of H40L is different. The mass of central object for 
H40L increases faster than any other models. This means that the accretion rate is highest 
among the models and leads to the shortest neutrino burst. Since the mass becomes large 
soon after bounce, the density and temperature inside the proto-neutron star become high 
and they cause the quick rise of average energies and luminosities as discussed above. This 
behavi or is in accord with the a nalyses for their numerical results of the collapse of 4OM0 



star bv ILiebendorfer et all (120041 ) 
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4.4. Profiles of density and temperature 

We examine the profiles of three models during the evolution to discuss the origin of 
differences in the black hole formation and the neutrino signals. 

In order to clarify the origin of different histories of accretion rate, we make comparisons 
of density profiles for W40L, T50L and H40L at the time when the central density reaches 
10^^ g/cm^ and at 30 ms after bounce in Fig. M In the upper panel of Fig. [9l we see the 
same profiles at center (M5 < IMq) when we compare at the same central density during 
the collapse, though we have seen the differences among the initial models as a whole in Fig. 
[H This similarity of central cores implies the similar dynamics around core bounce. In fact, 
the density profiles of central part (M5 < I.2M0) at 30 ms after bounce are quite similar as 
seen in the lower panel of Fig. [91 We note that profiles of temperature and electron fraction 
are also similar in the central part. 

It is remarkable to see the difference of densities in outer part (M5 > IM©) among 
the three models. Since the density at each point determines the time scale of free fall 

\jyJ~Gp) onto the central object, the accretion rate (M = ArKV^pv) is proportional to 
^ p^/^, therefore, the density difference is the origin of different accretion rates. In Fig. O 
the higher density in the outer profile of H40L informs that the accretion rate is higher up to 
the point at M5^2.4M0, where there is a crossing with the two other models. The profiles 
of T50L and W40L are rather similar each other though there are slight differences and a 
crossing at M5^2.1M0. 

The difference in the densities of outer layers remains during the evolution of proto- 
neutron star up to the black hole formation and leads to the different curves for proto-neutron 
star masses as we have seen in Fig. [8l At 30 ms after bounce, proto-neutron stars having 
~1.4M0 (in baryon mass) are formed and the rest of outer layers are accreting onto them in 
the free fall time. As they evolve, the baryon mass of proto-neutron star of H40L becomes 
larger than those of W40L and T50L because of the higher accretion rate. The cases of 
W40L and T50L in Fig. [8] are almost the same each other because of the similarity of density 
profiles. In the case of Shen-EOS, the cases of W40S and T50S in Fig. [8] behave similarly 
each other for the same reason. The increase in T50S at ^1 s becomes slightly slower than 
that in W40S refiecting a crossing in the density profiles similar to the one discussed above. 

In Figure [ini we compare the profiles of density and temperature at 300 ms after bounce 
for W40L, T50L and H40L. The central density for H40L is higher than those for W40L and 
T50L due to the larger baryon mass in proto-neutron star for H40L as seen in Fig. [8l The 
corresponding temperature for H40L is higher than others because of a more compact profile 
of proto-neutron star. We note that a kink at ^100 km in the density profile of H40L 
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corresponds to the shock position in accreting matter. The profiles of W40L and T50L are 
similar each other, though the profile in T50L is slightly more compact having higher density 
and temperature than those in W40L. 

The compact profile in H40L leads to the higher average energies and luminosities in 
neutrino emission than those in the two other models. The higher average energy of z/^/^ 
in H40L arises from the higher temperature in H40L than others. The higher luminosities 
in H40L is also caused by the higher temperature and the compactness of central object. 
Since the accretion luminosities of and are determined by the liberated energy in the 
gravitational potential of central proto-neutron star, a more massive and compact object 
leads to larger luminosities. In contrast, the neutrino energies and luminosities of W40L 
and T50L are similar each other because of the similar history of accretion and resulting 
evolution of proto-neutron stars. The same argument applies for the comparisons between 
W40S and T50S, in which the average energies are very close. The luminosities of T50S are 
somewhat smaller corresponding to the lower accretion rate at the late phase ^0.6 s after 
bounce. 

To summarize, the density profile in progenitors determines the accretion rate in black 
hole forming collapse and aflFects the evolution of proto-neutron stars toward the black hole 
formation. The difference in proto-neutron star profiles in turn leads to the diflFerent features 
in neutrino emission. The characteristics of neutrino emission refiect the minute differences 
in density profiles of progenitors while the infiuence from the equation of state is more 
significant. 

4.5. Differences in neutrino signals 

Here we discuss the characteristics of neutrino bursts in the five calculated models aiming 
at discrimination of the effects from progenitors and EOSs. We summarize the numerical 
results regarding the black hole formation in Table [H 

In order to probe the EOS, the time duration of neutrino burst is decisive. If the 
accretion rate of the progenitor can be determined, the measurement of the duration enables 
us to infer the critical mass of proto-neutron stars and to extract the stiffness of the EOS. 
In the current study, Shen-EOS is stiff and provides large critical masses of ^2.7M0 in 
baryon mass (^2.3M0 in gravitational mass, see Tabled]) for accreting proto-neutron stars. 
This large mass corresponds to the duration till the black hole formation over 1.3 s and is 
significantly longer than the cases of Lattimer-Swesty EOS as seen in Fig. [81 Since Lattimer- 
Swesty EOS is soft, the critical mass is ^2.1M0 in baryon mass (~2.OM0 in gravitational 
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mass). The duration is ~0.5 s with weak dependence on the progenitors among T50L, H40L 
and W40L. If one can get the information of the time duration accurately, one can infer 
the critical mass through the calculated curves of mass increase. For example, if we assume 
that the time duration of neutrino burst is 0.3 s, the critical mass is in the range between 
I.8M0 and 2.IM0 estimated from the range between W40L and H40L in Fig. [51 Therefore, 
one can estimate the critical mass with an error of ^O.SMq from the time duration within 
the uncertainty in the density profile of the progenitors. It is to be noted that the critical 
masses for proto-neutron stars, which contain abundant neutrinos at finite temperature, are 
diflFerent from the maximum masses for neutrino-less and cold neutron stars. 

To reveal the differences due to progenitors, it is necessary to examine the detailed 
information of average energies and luminosities. When Shen-EOS is adopted, one has to 
wait for ~1 s to have enough diflFerences between W40S and T50S. In case of Lattimer- 
Swesty EOS, it is rather difficult to distinguish the differences between the progenitors since 
the duration of burst is too short. The neutrino signals in W40L and T50L are too similar 
each other to discriminate. H40L is different from W40L and T50L, on the other hand, and 
may be distinguishable if the data on the average energy of z/^/^ and/or accurate luminosities 
are available. 

In Figure dU we compare the energy spectra of emitted neutrinos at 300 ms after bounce 
for H40L and W40L. The shapes of energy spectra for and are similar between H40L 
and W40L, however, the peak values are higher for H40L. The peak position for z/^/^ for 
H40L is shifted for a higher energy, though the peak value is close to that for W40L. These 
differences appear as higher luminosities of all species and higher average energy of v^j^ (but 
similar average energies for and for H40L as seen at 300 ms in Fig {71 We compared 
the energy spectra at 300 ms after bounce also for T50L and W40L (not shown here). The 
energy spectra are similar each other except for slightly higher peak values for and for 
T50L. 

We remark here that the discussions on neutrino signals so far do not take into account 
the neutrino oscillations. One should take into account the changes of spectra during the 
propagation of neutrinos in outer layers. The observational aspects of the current numerical 
results considering the effects of neutrino oscillations are now under investigation and will 
be published elsewhere. 
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5. Discussions 

In this section, we discuss the perspective of the current work regarding the progenitor 
and the equation of state. We also remark on the observational chances of neutrino burst 
and implications to astrophysical phenomena with the black hole. 

The present study on the black hole formation and the associated neutrino burst relies 
on the presupernova configurations of massive stars. Further studies with updated models 
of progenitors are definitely warranted. Necessary condition to have this type of phenomena 
is the appropriate density distribution to cause the intense accretion. In the three models 
(T50, H40 and W40) adopted in the current study, the size of iron core is large (~2M0) and 
the dense region ex t ends u p to M5 ^SMq. This is diflFerent from the case of a I5M0 star by 



Woosley fc Weaverl ( 119951 ). in which the iron core is smaller and the density slope is steep 
already at M5 ^1.5M0, for example. It is necessary to study further the density profiles of 
various progenitors and its relation to the accretion in order to examine whether the size of 
iron core is essential. It may be enough to have the fiat profile of density so that the density 
remains high ('--^lO^ g/cm^) at M5=2-3M0, which corresponds to the critical mass of proto- 
neutron stars. If the density profile is steep and the density is too low at the corresponding 
region in 4O-5OM0 stars, the black hole formation may take much longer (^10-100 s) than 
the current case. Such events with longer neutrino bursts may look similar to the ordinary 
supernovae of ^15M0 stars in terms of neutrino signals. 

Regarding the uncertainty of progenitors, o ne should pay atten ti on to the mass loss as 



well a s the rotation. The progenitor models in IWoosley fc Weaverl (119951 ) and [Hashimoto 



( I1995I ) are obtained by the presupernova evolution without the mass loss. The inclusion of 
the mass loss leads to the smaller mass at the precollapse stage even from the large initial 
mass (ex. 5OM0) at the main sequence. The iron core mass and the density distribution may 
end up to be similar to those in the progenitors of about 2OM0. In this regard, metal poor 
massive stars are more likely to be the progenitors for the black hole formation as studied in 
the current work. In fact, T50 adopted in the current study is the one with zero metallicity 
that forbids the mass loss. Its large central core results in the necessary intense accretion, 
which in turn leads to the quick black hole formation with short neutrino bursts. Whether 
the change of evolution toward the solar metallicity may cause the suppression of the current 
scenario remains to be seen along with the systematic studies on the metallicity dependence 
of progenitors and the careful studies on the mass loss. 

We study the black hole formation under the spherical symmetry. This geometrical 
assumption limits ourselves to the non-rotating massive stars among varieties in rotation 
rates. One should note, however, that the spherical symmetry is a good approximation if 
the rotation is not significantly fast. It is challenging to study the collapse of rapidly rotating 
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massive stars, which are possible progenitors for hypernovae and gamma ray bursts. One 
has to perform the fuUy general relativistic calculations in multi-dimensions to follow the 
evolution toward the black hole formation. Such challenges are made in the numerical simu- 
lations of general relativistic hydrodynamic s but with simple treatmen t of neutrino transfer, 
microphysi cs (ex. EOS) and progenitors (jShibata fc Sekiguchil l2005l : ISekiguchi fc Shibata 



20051 . 



20071 ). Multi-dimensional simulations with proper treatment of neutrino transfer are 



awaited for to predict the features of neutrino emission. 

The rotating collapse to the black hole and the following evolution is also interesting from 
the aspect of nucleosynthesis. In the context of coUapsar scenario, the ejection of neutron- 
rich material in the jet from the surroundings of black hole and the neutrino-driven winds 
from the accretion disk around the black hole are prop osed to be the possible site for the 



heav y elements production including r-process elements ( jSurman et al.ll2006l : iFujimoto et al. 



20071 ). Although we performed the numerical simulations in the limit of no rotation, one 
may use the current results to guess the time duration till the black hole formation and 
the accretion rate of material, which are influential in the scenario of nucleosynthesis. The 
time duration obtained in the current work is a lower limit since the maximum mass of 



rotating proto-neutron stars is larger than the non-rotating ones (ISumiyoshi et al.lll999l ) and 
the rotation supports accreting material away due to the centrifugal force. The accretion 
rate in the current study is an upper limit, accordingly. Large accretion rates in the disk 
formed around black hole, which are assumed in the studies of nucleosynthesis of r-process 
(ISurman et al.l l2006l ) , may be limited from the density proflle of progenitors as we have 
pointed out. 

The observational possibility to detect the neutrino bursts toward the black hole for- 
mation from massive stars deserves to be mentioned. Occurrence rate of the phenomena 
(non-rotating black hole formation and its short neutrino burst) as in the current scenario 
is proportional to the fraction of massive stars having large masses for intense accretion 
without signiflcant rotation. This depends on the initial mass function and rotation rates 
with large uncertainties, but such stars comprise a certain fraction (^30%) of massive stars, 
which can be estimated from the Salpeter's initial mass f unction with an a s sump tion on the 
mass threshold (^ 25M0) for the black hole formation. iKochanek et al.l (120081 ) estimated 
from the observations done so far that the rate of black hole formation could be comparable 
to the rate of normal core-collapse supernovae. Observational analyses of explosion energies 
and Ni productions for various supernovae indicate that there have been observed at least 
several faint supernovae (SN 1997D and SN 1999br, for example) from massive stars beyond 
20Mr:^ without signif icant rotations, which may be akin to the event considered in this paper 
jNomoto et al.ll2007l ). 
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The planned surv ey of disappearance of massive stars within a distance of 10 Mpc 
( Kochanek et al. 20081 ) will find both dim supernovae and collapses with no optical display 
that lead to the black hole formation. If they occur in the Galaxy or nearby galaxies, the 
information on progenitor will be obtained by the optical records before it disappeared. In 
addition, thousands of neutrinos will be detected at the terrestrial facilities. Neutrino data 
combined with optical ones will provide us with the information on the dynamics of black 
hole formation and constrain the properties of dense matter and the progenitor. In the search 
of supernova neutrinos by the Super-Kamiokande detector, neither neutrinos associated with 
black hole formation n or ordinary supern ova neutrinos are detected during the period of their 
survey unfortunately (jlkeda et al.l 120071 ). 



The numerical result on the massive star with zero metallicity (T50) suggests that the 
black ho le formation in t h e current sc enario may occur for Pop III stars in the early Universe 
(see also iNakazato et al.l (120061 . 120071 )). Since the massive stars become more popular in the 
metal poor epoch, the neutrinos from black hole formation in the early stage of galaxies 
may contribute to the relic signal as back ground neutrinos. Although the signal is short, 
higher energies and luminosities may help to enhance this contribution and they may not 
be negligible. This argument again depends on the fraction of massive stars resulting the 
current scenario through the evolution of galaxies. 

Before we close the discussion, we comment on the choice of equation of state. Although 
density and temperature become extremely high enough to have new degrees of freedom 
such as hyperons and quarks in the evolution toward the black hole formation, we adopt the 
two sets of EOS within the nucleonic degree of freedom in the current study to assess the 
progenitor dependence. Studies on the modifications of the evolution of black hole formation 
and the neutrino signal due to hyperon mixtures an d the appearance of qu ark phase are now 
under way and will be reported in separate papers (INakazato et al.ll2008l ). 



6. Summary 



We study the black hole formation and the associated neutrino emission from the gravi- 
tational collapse of massive stars to clarify the depen dence on the progenito r models following 
the basic study of the phenomena in the first paper (ISumiyoshi et al.ll2007l ) of the series. We 
perform the sets of numerical simulations of general relativistic z/-radiation hydrodynamics 
under the spherical symmetry to follow the collapse of massive stars of 4OM0 and 5OM0 
adopting different models of the presupernova evolution. We clarify that the characteristics 
of evolution of proto-neutron stars toward the black hole formation and the associated short 
burst of neutrinos are found similar in the numerical simulations for the three progenitor 
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models. 

We find that the different progenitor models may result in different rates of matter 
accretion onto the proto-neutron star born in the collapse of massive stars. The different 
accretion rates can lead to different histories of increasing mass of proto-neutron stars toward 
the critical mass and to different durations till the black hole formation. In order to extract 
the properties of equation of state from the neutrino burst toward the black hole forma tion 



from massive stars, as proposed in the previous studies (ISumiyoshi et al.l l2006l . 120071 ). it 
is necessary to constrain the uncertainty from progenitor models. In the comparison of 
the three progenitors we studied, the change of time duration is found not so large as the 
difference originating from the two sets of equation of state. The change of time profiles of 
average energies and luminosities of neutrinos due to progenitors is found to be minor as 
well. Note, however, that the present study is a first trial exploration and more systematic 
investigations are obviously necessary to see if our conclusion is generic or not. 

We demonstrate that there is a relation between the density profile of progenitor and the 
accretion history to cause the difference in black hole formation. A higher density at M5=2.0- 
3.OM0 in a progenitor results in a shorter free fall time and a higher accretion rate. A flat 
density profile of core in massive stars is a key factor to cause the rapid formation of black 
hole formation in ^1 s together with the characteristic neutrino burst. When the density 
gradients are different as in the comparison between W40L and H40L, the density profiles 
may be distinguished from the detailed information of neutrino detections. This relation of 
the density profile of progenitors to the accretion rates is helpful to discuss the evolution 
after the collapse of massive stars in other stellar models and may be useful to infer the 
accretion dynamics of rotating collapse to the black hole formation and associated explosive 
phenomena such as coUapsars. Further studies on rotating collapse of massive stars with 
intense accretion are called for to clarify whether the outcome may lead to possible success 
of modeling gamma ray burst and/or nucleosynthesis of heavy elements. 
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Table 1: Summary of calculated models 



Model Progenitor Mprgg Mpe EOS Mg^"^ [Mq] M^"^ [Mq] t^H [s; 



W40S 


WW95^ 


40 


1.98 


Shen 


2.66 


2.38 


1.35 


W40L 


WW95 


40 


1.98 


LS 


2.10 


1.99 


0.57 


T50S 


TUNO?*^ 


50 


1.88 


Shen 


2.65 


2.33 


1.51 


T50L 


TUN07 


50 


1.88 


LS 


2.11 


2.01 


0.51 


H40L 


H95^ 


40 


1.88 


LS 


2.17 


2.08 


0.36 



^WW95: IWQQslev &: Weaverl (|l995h 

^TUN Q7: lUmeda fc Nomotol (|2QQ5l ): IXominaga et"al] (|2QQ7h 
^H95: iHashimotol (|l995h 



Note. — Mprog and Mp^ is the mass of progenitor and the iron core mass, respectively. M^'^^ and 
y^max ^^le masscs of proto-neutron star just prior to the re-cohapse in baryon and gravitational masses, 
respectively. tsH is the time at the black hole formation measured from the core bounce. See the main text 
for details. 
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Fig. 1. — Density, electron fraction and temperature distributions of progenitor models are 
shown as a function of baryon mass coordinate by solid, dashed and dash-dotted lines for 
W40, T50 and H40, respectively 
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Fig. 2. — Radial trajectories of mass elements of the core of SOMq star as a function of time 
after bounce in model T50S. The location of shock wave is displayed by a thick dashed line. 




tiim[sw] 



Fig. 3. — Radial trajectories of mass elements of the core of 5OM0 star as a function of time 
after bounce in model T50L. The location of shock wave is displayed by a thick dashed line. 
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Fig. 4. — Average energies and luminosities of (solid), (dashed) and z/^/^ (dash-dotted) 
for model T50S are shown as a function of time after bounce. The results for model W40S 
are shown by thin lines. 
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Fig. 5. — Average energies and luminosities of (solid), (dashed) and v^j^ (dash-dotted) 
for model T50L are shown as a function of time after bounce. The results for model W40L 
are shown by thin lines. 
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Fig. 6. — Radial trajectories of mass elements of the core of 4OM0 star as a function of time 
after bounce in model H40L. The location of shock wave is displayed by a thick dashed line. 
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Fig. 7. — Average energies and luminosities of (solid), (dashed) and v^j^ (dash-dotted) 
for model H40L are shown as a function of time after bounce. The results for model W40L 
are shown by thin lines. 
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Fig. 8. — Enclosed baryon mass in the central object (proto-neutron star) as a function 
of time after bounce for models W40S (long-dashed), W40L (sohd), T50S (dotted), T50L 
(dashed) and H40L (dash-dotted). 
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Fig. 9. — Density distributions for W40L, T50L and H40L are shown as a function of baryon 
mass coordinate by solid, dashed and dash-dotted lines, respectively. Upper and lower panels 
show the profiles at the time when the central density reaches 10^^ g/cm^ and at 30 ms after 
bounce, respectively. 
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Fig. 10. — Density (upper panel) and temperature (lower panel) distributions for W40L, 
T50L and H40L at 300 ms after bounce are shown as a function of radius by solid, dashed 
and dash-dotted lines, respectively. 
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Fig. 11. — Energy spectra at 300 ms after bounce are shown as a function of neutrino energy 
for i^e, and v^^jr by solid, dashed and dash-dotted hues, respectively, for model H40L 
(thick) and W40L (thin). 



